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ABSTRACT 

Context. Microlensing is a unique method to probe low mass exoplanets beyond the snow line. However, the scientific potential of 
the new microlensing planet discovery is often unfulfilled due to lack of knowledge of the properties of the lens and source stars. The 
discovery light curve of the super Earth MOA-2007-BLG-192Lb suffers from significant degeneracies that limit what can be inferred 
about its physical properties. 

Aims. High resolution adaptive optics images allow us to solve this problem by resolving the microlensing target from all unrelated 
background stars, yielding the unique determination of magnified source and lens fluxes. This estimation permits the solution of our 
microlens model for the mass of the planet and its host and their physical projected separation. 

Methods. We observed the microlensing event MOA-2007-BLG-192 at high angular resolution in JHKs with the NACO adaptive 
optics system on the VLT while the object was still amplified by a factor 1.23 and then at baseline 18 months later. We analyzed and 
calibrated the NACO photometry in the standard 2MASS system in order to accurately constrain the source and the lens star fluxes. 
Results. We detect light from the host star of MOA-2007-BLG-192Lb , which significantly reduces the uncertainties in its char- 
acteristics as compared to earlier analyses. We find that MOA-2007-BLG-192L is most likely a very low mass late type M-dwarf 
(0.084+°;°^ M Q ) at a distance of 660+^° pc orbited by a 3.2^; 2 M e super-Earth at 0.66+°;^ AU. We then discuss the properties of this 
cold planetary system. 

Key words, -techniques: microlensing, image processing -exoplanets: individual: MOA-2007-BLG-192Lb, -stars: low-mass stars, 
late-type 



1. Introduction 

Gravitational microlensing provides a unique window on extra- 
solar planetary systems with sen sitivity to cool planets, part icu- 
larly those of low mass ([Bennett & Rhie 1996; Beaul ieu et al. 
200 6]|Gould etal |2006[ |Bennett et al.|2008t|Kubas et al.|2 008 ; 



Sum i et al.| |2010) that are currently well beyond the reach of 
other methods. Microlensing is also sensitive to planets orbiting 
very faint stars and hence spectral types not routinely examined 
with other techniques. In general it is a powerful tool to study the 
Galactic planetary population as a whole ( [Cassan et al.|2012| ). 



* Based on observations under ESO Prog.IDs: 279.C-5044(A) and 
383-C.0495(A) 



Microlensing occurs when a foreground (lens) star passes 
close to the line of sight towards a background (source) star. The 
gravity of the foreground star acts as a magnifying lens, increas- 
ing the apparent brightness of the background star as it gets close 
to the line of sight. A planetary companion to the lens star will 
induce a perturbation to the microlensing light curve with a du- 
ration that scales with the square root of the planet mass, lasting 
typically a few hours for an Earth to a few days for a Jupite r 
( [Gould & Loebl[T992l |Mao & Paczynsk!l[T99Tj |Liebes|[T^64l ). 



Hence planets are now routinely discovered by dense photomet- 
ric sampling of ongoing microlensing events. The inverse prob- 
lem, finding the properties of the lensing system from an ob- 
served light curve, is a complex non-linear one within a wide 
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parameter space. The planet/star mass ratio and projected star- 
planet separation can usually be measured with high precision. 
However in the absence of higher order effects such as parallax 
motion and/or extended source effects, in general there are no 
direct constraints on the physical masses and orbits of the plan- 
etary system. In the least information case, model distributions 
for the spatial mass density of the Milky Way, the velocity distri- 
bution of potential lens and source stars, and the mass function 
of the lens stars are used in a Bayesian analysis to derive proba- 
bility distributions for the masses of the planet and the lens star 
and their distance, as well as the orbital radius and period of the 
planet. 

With complementary high angular resolution observations, 
currently done either by HST or with adaptive optics, it is possi- 
ble to get additional strong constraints on the system parameters 
and determine masses to about 10%. This can be done by di- 
rectly measuring the light coming from the lens and measuring 
the lens and source relative proper motion (Benne ttet al.|2 006, 



2007] [20T0l |Gaudi et al. 12008} |Dong et al.|2009[|Janczak et al. 



2010) 



An extrasolar planet with a best- fit mass ratio of q ~ 2 x 10" 4 
was discovered in the microlensing event MOA 2007-BLG-192 
( [Bennett et al.||2008| ) found by the MOA collaboration toward 



the Galactic bulge, (J2000: RA, Dec ) = (18:08:03.8,-27:09:00). 
The best fit microlensing model shows both microlensing paral- 
lax and finite source effects. Combining these, we obtained the 
lens masses of Mi = 0.06+ q^M q for the primary and 3.3^M© 
for the planet. The incomplete light curve coverage of the plane- 
tary anomaly led to a significant degeneracy in the lens models, 
and the lack of strong constraints on the source size to a poorly 
determined Einstein radius. Together this resulted in rather large 
uncertainties in the physical parameter estimates of the system. 

Additional constraints are required to exclude competing mi- 
crolens solutions and to refine our knowledge of the physical 
parameters of the system. It is possible to constrain masses and 
parameters of the system thanks to high angular resolution imag- 
ing. Most microlensing events only provide a single parameter, 
the Einstein ring crossing time ts 9 that depends on the mass of 
the lens system M L , its distance D L , the source distance D s and 
their relative velocity. However, when the relative lens-source 
proper motion yu re i can be determined this yields the angular 
Einstein ring radius 6e = //reite- Moreover E is linked to the 
lens system mass by 



coming from the source and from the lens star to refine estimates 
of the parameters of the system. 

2. The data set 

We obtained JHKs measurements using the NACO AO system 

|2003| ) mounted on Yepun 



( [Lenzen et al. 2003 ; Rous se t et al. 



during the night 6/7 Sept. 2007, while the source star was still 
magnified by a factor of 1.23. AO corrections were performed 
on a natural guide StarBand observations with the S27 objective 
(27" x 27" FOV, pixelscale=0.02715") were conducted in jitter 
mode with multiple exposures at random offsets within 10" of 
the target. In the absence of suitable "empty" sky patch close to 
the target, this strategy was chosen to ensure an accurate esti- 
mation of the sky background and to filter out bad pixels. The 
second epoch(s) were obtained with the same observing strategy 
more than 22 months later with the event being at baseline, i.e. 
when the source was not magnified anymore. An overview of the 
NACO data set is given in Table [T] 

To perform absolute calibration of the NACO images we ob- 
tained 90 x 10 s dithered images in JHKs of the MOA-2007- 
BLG-192 field with t he Sirius simultaneous 3-band camera 
(Nagayama et al.|2003|) on the Japa nese/South African IRSF 1.4 
m telescope ( [Glass & Nagata|2000| ) at SAAO (non AO, ~ 8' x 8' 
FOV, pixelscale=0.45") on 29th of Aug. 2008, i.e. at a time when 
the event was at baseline. 



2.1. Reduction 

Following a "lucky imaging" approach we visually inspect each 
of the NACO raw images and remove the ones for which the AO 
correction was obviously poor. The remaining raw frames are 
then dark- subtracted with darks of exposure times matching the 
science frames, flatfielded with skyflats, median co-added and 
sky-subtracted using recipes from the Jitter/Eclipse infrared data 
reduction package by Devillardj ( |1997| [1999 ). To avoid border 
effects, we keep only the intersection of the different dither po- 
sitions of the co-added frames for our photometric analysis. 
The IRSF data, which was taken to gauge our NACO data, 
has been dark-subtracted, flat-fielded and sky subtracted using 
the on-the-mountain pipeline package for the SIRIUS camera 
( Nagayama et al.|2003| ). 



4G E D S 
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Therefore, since the distance of the source Dg is known from its 
magnitude and colors, Eq. ([I]) is a mass-distance relation for the 
lens star. Another constraint is needed to obtain a complete solu- 
tion to the microlensing event. This can be achieved by directly 
detecting light from the planetary host star (the lens). Combining 
this measurement with Eq. ([T]) and a mass luminosity relation 
will yield the mass of the lens. This has been done already for 
several microlensing eve nts where the system is composed of a 
star and a gaseous planet ([Bond et al.|2004[|Bennett et al.]2006 ; 



Gaudi et al. 2008 ; Janczak 



[ Udalski et al.p005[ |Dong et al.| 2009 
let al.|2010) . 

We observed MOA-2007-BLG-192 in JHK using adaptive op- 
tics on the VLT while it was still amplified by a factor of 1.23 
and again when the microlensing was over. Here, we combine 
the NACO JHK flux measurements at these 2 epochs with the 
color estimate of the source star ([Gould et al.|2010) and the mi 



3. Photometric Analysis 

As in our previous analysis of planetary microlensing event 
MOA-2008-BLG-310Lb ( [Janczak et aTpJIO) , we extract the 
photometry of NACO images using Starfinder ( [Diolaiti et al. 
2000 ). This tool is tailor suited to perform photometry of AO im- 
ages of crowded fields. It creates a numerical PSF template from 
chosen stars within the frame, which is then used for PSF-fitting 
of all stars in the field. To build our PSF reference we chose the 
star marked as "1" in Fig. [T] based on the following criteria. It 
is close to the target (within less than 4"), sufficiently bright but 
well within the linearity regime of the detector and common to 
all final reduced JHKs images of both epochs. Fig. [6] shows the 
JHKs images centered on the target for the first epoch of NACO 



crolensing model (Bennett et al. 



1 The LGSF, which in theory should have yielded better performance, 
was not available at that time. However the pro and contra of LGS vs 
NGS for us have to be evaluated on a case by case base, since in the 
crowded field of microlensing targets one often finds suitable NGS ref- 
erences which may give even better corrections than the LGSF accord- 
2008| ) to disentangle the flux ing to the ETC observation preparation software. 
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Fig. 2. Left: The (J-Ks, J) CMD in the 2MASS system of the MOA-2007-BLG-192 field combining the data from the IRSF 
(within 3' of target, black points) and NACO (within 18", blue points). In red the photometry of the measured lens+source flux 
at magnification A=1.23 is displayed together with the inferred decomposed fluxes of the source (green) and the lens (planetary 
host star, black). Overplotted are Marigo et al (2008) solar metallicity isochrones of ages log(Gyrs) = 9.00, 9.88, 10.15 at distance 
modulus of dm= 14.38 and estimated extinction of Aj = 0.72, Ak s = 0.29. Right: Same as above but for (H-Ks, H). 



Table 1. Log of JHKs NACO data. According to the Paranal 
night logs the Epoch 1 night was classified as photometric, 
whereas the Epoch 2 observations were taken in clear sky con- 
dition. We give the exposure time, modified Julian Date, airmass 
and measured full width at half max on the coadded frames. 



NACO Ks 



IRSF Ks 



Band 


n x Exp [s] 


MJD 


Airmass 


FWHM ["] 


Epoch 1 


J 


6x60 


54350.00781250 


1.005 


0.14 


H 


20x25 


54350.02734375 


1.023 


0.19 


Ks 


10x25 


54349.98828125 


1.002 


0.09 


Epoch 2 


J 


23x60 


55036.08593750 


1.015 


0.34 


H 


22x30 


55036.06640625 


1.034 


0.29 


Ks 


24x30 


55015.10156250 


1.088 


0.10 



and the PSF subtracted re siduals. The IRSF pho tometry catalog 
was created with DoPhot ( Schechter et al.|1993| ) 



3.1. Building a calibration ladder 

In order to build the calibration ladder, we use three data sets: the 
2MASS catalogue, photometry obtained at the IRSF telescope 
and then NACO data. It is necessary to use the intermediate step 
of IRSF observations because we have too few stars in common 
between 2MASS and NACO. Suc h ladder has been u sed already 
and described in the appendix of ( Janczak et al.|2010) . The three 
JHKs color systems are very close. For example, the color term 
for J-Ks colors between 2MASS and IRSF is only 0.01. Accurate 
calibration between 2MAS S and IRSF has been given in Kato 
|et al.| ( |2007T ). |Janczak et al.H2010) did not detect color terms be- 
tween NACO and IRSF filters. Therefore, our calibration ladder 
must only determine the zeropoint offsets. 
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Fig. 1. Left: Extract of IRSF Ks band image of MOA-2007- 
BLG-192 used to calibrate the NACO photometry of the 18" x 
18" large intersection fov of the coadded NACO frames in Ks 
band (right). MOA-2007-BLG-192 is marked with the half cross 
hair. The stars annotated with "1" and "2" serve as PSF-reference 
and photometric zeropoint calibrators. Furthermore these two 
stars are common to all bands and epochs. The bright stars north 
of the two references are either too crowded, in the non-linear 
regime or too far away from the target. 



We first perform the astrometry of the IRSF images with 
respect to the online 2MASS catalog using GAI A/Sky cat and 
WCSTools. Then, using only stars marked as AAA (highest 
2MASS quality flag) in the JHKs bands we crossmatch the com- 
mon stars to compute the photometric transformation between 
the two catalogs by sigma clipping, demanding an astrometric 
accuracy of the match of better than 0.6". To minimize the ef- 
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Fig. 3. Same as Fig.^but for (J-H,J). 



Table 2. JHKs NACO photometry for MOA-2007-BLG-192, i.e. 
lens + source (no dereddening applied). The absolute photometry 
error budget is composed by adding in quadrature the errors on 
the zeropoint, the formal error reported by Starfinder and the 
background error as estimated from the scatter between epoch 1 
and epoch 2 comparison stars. For epoch 1, 7, and H bands, we 
adopt the background error estimate as derived from the K band, 
since the poor epoch 2 quality in J and H would overestimate the 
epoch 1 errors. 



Band / H 




NACO Epoch 1 


calibrated against IRSF 


19.209 + 0.043 18.281 +0.042 


17.948 + 0.035 


NACO Epoch 2 


calibrated against IRSF 


19.324 + 0.073 18.548 + 0.112 


17.989 + 0.038 


A Epochs 0.115 + 0.085 0.267 + 0.120 


0.041 + 0.052 


aligned with respect to Epoch 1 


19.283 + 0.071 18.498 + 0.087 


18.011+0.042 


A Epochs 0.074 + 0.083 0.217 + 0.097 


0.063 + 0.055 



Aj = 0.72 + 0.10, A H = 0.46 + 0.10, A K = 0.29 + 0.10, which for 



this line of sight is consistent with extinction maps from Schlegel 
|eTaLl ( [T998l ) and |Marshall etaL] ( [2006] ) 



feet of source confusion and blending contamination we cut off 
at magnitude 13 for the 2MASS reference stars and sum up the 
flux of close neighbors for the IRSF sources to account for the 
much coarser pixel scale of the 2MASS catalog. 

The PSF reference star is contained in the IRSF catalog, as 
well as star "2" (Fig. [I]). We examine their long term photomet- 
ric stability in the OGLE database and find that over more than 
seven years both stars are stable (in the optical /-band) at levels 
of < 1%, which makes them well suited as zeropoint calibrators 
of our NACO field. While we adopt star "1" as the primary pho- 
tometric calibrator since star "2" is more crowded, we determine 
zeropoints from both stars as a consistency check. To account for 
the different plate scales between NACO and IRSF we sum up 
the flux of all the NACO sources which are contained within the 
IRSF PSF. We note that observing conditions (sky transparency 
and atmospheric coherence times) for the second epoch data set 
were inferior to the epoch 1 measurements and the uncertainties 
in the absolute zeropoints of epoch 2 are therefore larger. Since 
we are mainly interested however in the relative photometry of 
the two epochs we can align the epoch 2 photometry with respect 
to more accurately calibrated epoch 1. Table |A| summarizes this 
way of determining the transformations to calibrate the NACO 
data with respect to the 2MASS system and Table [2] shows our 
derived photometry for MOA-2007-BLG-192. 



4. Results 

In Fig. [2] and [3] we present the color-magnitude diagrams for 
the combined IRSF and NACO (epoch 1) data. To estimate the 
interstellar extinction, we first determine the position of the red 
clump center by taking the median of the distributions in color 
and magnitude inside a window centered on a first guess esti- 
mated position. Then we fit the tip of the Red Gi ant Branch as 
given by the isochrones of Mari go et al.| ( |2008| ) adopting the 
distance modulus dm = 14.38 + 0.07 as found for the MOA- 
2007-BLG-192 field by |Bennett eFaL] ( [2008] ). With a best fit 
age of log(Gr) = 9.88 we find for the extinction coefficients: 



4.1. The case for a luminous lens I: NACO-only 

The standard general microlens light curve model is given as 

F(t) = F s A(t) + F B , (2) 

where F is the measured flux at the telescope, F$ is the intrinsic 
unmagnified source flux, A(t) the time dependent magnification 
given by the lens model and the blend flux F B = F L + ^Background* 
which contains the lens flux F L and ^Background the flux of any 
unrelated field stars within the aperture, FBackground- While the 
source flux F$ can be determined with high precision from the 
light curve modeling of the non-AO data, given a large magni- 
fication gradient, the background term normally dominates over 
the lens term in seeing-limited photometry of typically crowded 
Galactic Bulge fields of microlensing. Hence the benefits of high 
spatial resolution imaging are obvious. Reducing or eliminating 
the contribution of contaminating background sources provides 
a better estimate of the lens flux and so finally of the phys ical 
characteristics of the lens system. In Janczak et al. ( |2010| ) the 
lens flux could be estimated by comparing a single NACO AO 
epoch with an excellent seeing-limited light curve in the same 
passband from which the source flux had been previously de- 
termined with good accuracy. For MOA-2007-BLG-192Lb we 
have no such light curve in the NACO passbands but a well de- 
termined measurement of the source flux in the / band (Cousins 
system), which we can transform into the expected source flux 
for JHKs bands. Note that while in theory our two point NACO 
"light curve" can be used to solve Eq. ^ for the lens and source 
fluxes directly, the resulting uncertainties are very large (Dong 
|et al.|2009] ) due to very small magnification "lever arm" for our 
event and so following the path of Janczak et al. ( 2010| ) is much 
more accurate. 

First, however, the two epochs can be used as follows, with- 
out the knowledge of the source flux, to check whether there 
is an indication that light from the lens is detected. The ex- 
pected magnification gradient between the tw o NA CO epochs 
based on the best-fit model of Benn ett et al.| ( |2008j > is Am = 
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Fig. 4. Relative photometry of the NACO epochs for each band calibrated / aligned to Epoch 1 based on stars within 4" of MOA- 
2007-BLG-192. The black points are constant stars, so their scatter gives an estimate of the error. The red dashed line marks the 
expected magnitude difference of the target assuming no light from the lens is detected. The red diamonds show the flux changes of 
the target, which are inconsistent with such a scenario at a 3 - cr level for J and even with higher significance for our best data set 
in K band. The poor quality of the epoch 2 H band data as evident by the large scatter does not provide any useful constraint. Table 
[2] gives the derived values for all bands and different choices of the field of comparison stars. 



0.230 ±0.015 mag. Note that this gradient is basically the same 
for all competing planetary models, since the first epoch was 
taken close to the baseline of the event, where the single lens 
approximation describes the data very well. If the lens is dark 
and no unrelated source is contaminating our photometry (see 
Sec. |4.3| ) we then would expect to measure this difference in the 
relative photometry of the two epochs in each band. Since the 
quality of epoch 1 is superior we choose epoch 1 as reference 
to which we align epoch 2. We compare the photometry be- 
tween the two epochs for each band using 3 different alignment 
procedures. First we compare the derived absolute photometry 
(with r espe ct to 2MASS using the calibration ladder described 
in Sec. [37]}. 

Then we align epoch 2 with respect to (calibrated) epoch 1 
using all common stars within 4" (to minimize effect of PSF 
variations) of the target. The resulting magnitude differences for 
the target and the absolute photometry values are summarized in 
Table [5] The difference between the epochs is shown in Fig. [4j 
Regardless of the alignment method used, for all bands except 
H (the set with the poorest epoch 2 data quality), the measured 
difference is less than in the case of a dark lens, albeit with dif- 
ferent levels of significance. For K band, the best data set, a dark 
lens is inconsistent with the measurement at 2<x for the absolute 
alignment and at 3<x for the relative alignment. The results for 
/ band are also inconsistent with a dark lens, but in this case at 
slightly less than 2 cr. This motivates a more careful examination 
of the evidence for a luminous lens. 



4.2. The case for a luminous lens II: NACO+IRSF+Optical 

A more powerful test for the presence of "excess light" (in addi- 
tion to that of the source) is possible by combining NACO, IRSF, 
and optical data. To maximize sensitivity, we will work entirely 
with uncalibrated data. This will eliminate any errors associated 
with calibration relations, extinction estimates, and color-color 
relations. The remaining errors, which are either measurement 
errors or intrinsic scatter, are both small and easy to measure. 

We begin by constructing a color-color diagram that com- 
bines optical V and / data from OGLE-II with K data from IRSF 



(see Fig. [5]). The OGLE-II data are used because this is the sys- 
tem in which [Gould et al.| ( |201Q| ) measured the color and magni- 
tude of MOA-2007-BLG-192S, 

(V - I) s = 2.36 ± 0.03; /, = 21.45 ± 0.05. (3) 

The OGLE-II data are in fact calibrated, but that is incidental: 
the important thing is that the optical color is measured in this 
system. The K magnitudes are constructed directly from IRSF 
photometry fluxes K = 22.155 - 2.5 log(flux). The zero-point 
constant is chosen for convenience to be similar to the calibration 
constant, but this constant does not enter the calculation in any 
way. In particular, the data remain uncalibrated because there is 
no color term. 

The open black circles are all the astrometric matches that 
meet a strict crowding criterion. A color-color relation (red) is 
derived by fitting the points in the range 2.0 < (V - 1) < 2.4 with 
3 cr rejection (green filled points). This choice of interval will be 
justified below. The relation is: 

(I - K) = Zi + Z 2 [(V -I)- 2.36]; 
Z x = 2.757 ± 0.008; Z 2 = 0.943 ± 0.039 (4) 

with a scatter of 0.080 mag. 

The vertical lines represent the best fit and error bar of the 
Gould et al. (2010) optical-color measurement. From Eq.s ^ 
and ([4]), the best estimate of the source magnitude at the first 
epoch (when the source was magnified by A = 1.23) is therefore 

K s = [I s - 2.5 log(A)] - Zi - Z 2 [(V - /), - 2.36] = 18.468 (5) 

We discuss all the errors in this estimate below. 
We now compare this with the baseline flux as measured by 
NACO and transformed to the IRSF system: 

^baseJRSF = 

^compJRSF + 2.5 log^comp^NACo/^bascNACo) = 17.948 (6) 

where ^ C omp,iRSF = 14.261 ± 0.016 is the IRSF magnitude of the 
comparison star, and F comp ,NACo = 88066 ± 37 and F base,NAC0 = 
2952.2 ± 12.1 are the NACO fluxes of the comparison and base- 
line stars, respectively. 
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Fig. 5. Empirical VIK instrumental color-color relation. Lower 
panel: Open black symbols show all astrometric matches (that 
survived a strict crowding criterion) between OGLE-II V/I data 
and instrumental IRSF K data. The green filled points show 
those used in the fit (red line), which are restricted to 2.0 < 
(V - I) < 2.4 and exclude 3 cr outliers. Vertical magenta lines 
show MOA-2007-BLG-192S (V - I) color measurement by 
Gould et aL] ( |2010| ). Upper panel shows the VI CMD positions 
of stars used to determine the VI K relation. They are all giants 
in or near the clum p, whereas the source ( magenta point) is a 
dwarf. Nevertheless, |Bessell & Brett| ( [T988] ) show that the VI K 
relations are essentially identical for giants and dwarfs in this 
color range (see text). 



The baseline flux is clearly larger, K s - ^ase = 0.520, The 
question is, how large is the error in this difference? Eq. ([5]) has 5 
identifiable sources of error. First, the error in I s is 0.05, but the 
error in [I s - 2.5 log(A)] is actually smaller than this by a factor 
0.57 (see Eq. 10 of |Janczak et al.|20T0| ). Second, the error in Z\ 
is 0.008. Third the error in the final term is Z2 times the error in 
(V - I) s , i.e., 0.056. 

Fourth, we are using the VI K color-color relation of the field 
stars observed by IRSF as a proxy for the VIK color-color re- 
lation of the source. However, the source is a dwarf, while the 
field stars ar e all giants. Now, according to Fig. 1 of Bessell & 
Brett] ( |1988| ) these VIK relations of giants and dwarfs are virtu- 
ally identical for (V-K) < 3.0 [(V-I) < 1.3], and diverge only 
very slowly at redder colors. The dereddened color of the source 
is (V - 7)o = 1.24 + 0.06, so the entire 1 cr error range lies within 
the "same relation" region. And again, the relations diverge only 
very slowly at redder colors. We note, however, from Fig. 2 of 
Bessell & Brett (1988), that the divergence is extremely rapid in 



(J - K). This is the principal reason that we conduct this test in 
K rather than /. 

Finally, the VIK relation in Fig. [5] exhibits a scatter of 0.08 
mag. If this scatter is attributed to measurement errors, then the 
effect is very small. T here is some reason to expect that this is the 
case because Fig. 1 of Bessell & Brett| ( |l988 ) shows almost zero 
scatter. However, the bulge star population may be more diverse 
than the local one. For the moment we assume it is intrinsic, and 
that this scatter in the observed giant VIK relation also applies 
to dwarfs. Then, the error in K S (A = 1.23) is [0.028 2 + 0.008 2 + 
0.056 2 + 0.08 2 ] 1/2 = 0.102. 

There are three errors contributing to Z^base, which are listed 
below Eq. ([6]). Their sum in quadrature is 0.017. As seen in sec- 
tion 3 the JHKs bands from 2MASS, IRSF and NACO are very 
close, with negligible color terms in the transformations between 
the different systems. Next we note that photometry on AO im- 
ages is notoriously difficult due to gradients in the PSF. This 
effect is hard to quantify. However, Fig. [6] shows that it must be 
quite small. The left (right) column shows the image before (af- 
ter) PSF subtraction for the first epoch, for /, H, K, respectively. 
The K image in particular, looks extremely clean. We neverthe- 
less conservatively estimate a 0.03 photometry error due to PSF 
gradient. This yields a total error on ^ base of 0.034 mag. 

Finally, we note that even if the scatter in Fig. [5] were due to 
measurement noise, rather than intrinsic scatter (and so should 
not have been included in the error in K S (A = 1.23), it would 
then contribute to the error in ^ase through AT C omp,iRSF- Hence, 
the impact on the final error would have been identical. 

We therefore finally derive our estimate of the excess mag- 
nitude at baseline: 



K s - ^ba 



0.520 + 0.108 



which is a 4.8 cr detection of additional light. It is either a blend 
aligned to better than 0.1 arcsec with the source star of the mi- 
crolensing event or light coming from the lens star. 

4.3. Is the blended light from the lens star? 

The mean density of stars of comparable brightness and color 
±0.20 mag to the detected blend is less than 0.2 per arcsec 2 as 
derived from our best / sharpest data set, the Ks band of epoch 
1. Given the image quality of 0.09" FWHM, this conservatively 
implies a probability of less than 2% for the blend being unre- 
lated to the microlens event. Another possibility to consider is 
that the blend stems from a companion to the source star. Close 
companions with periods < 100 d can be ruled out by the xal- 
larap signal limits in the light curve and very wide separation 
companions > 700 AU would be resolved in the Ks NACO 
data. This still leaves a large range of allowed separations but 
taking into account the color difference the possible fraction of 
low mass secondaries should not be larger then 8% according to 
Duquenn oy & Mayor| ( |199l] ). However only future AO or HST 
images, when the source and lens will have moved sufficiently 
far apart to be spatially resolved, will be able to securely rule out 
such a scenario. 

4.4. Source star constraints 



To compare the (V - I)o = 1.24 ± 0.06 color found by |Gould 
|et al.| ( |2010| ) to the NIR bands of this study we transform this 
V — I color to / - K in the 2MASS system using first the dwarf 
color table of |Bessell & Brett| | [l958) to find (J-K) 0fS = 0.73 and 
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NACO Epoch 1 zoom (3.3" x 3.3") 
original after PSF subtraction 




Fig. 6. Left: Zoom (3.3" x 3.3" ) of Epoch 1 NACO images cen- 
tered on the target in H (top), / (middle), and K (bottom). Right: 
Images after PSF subtraction. The PSF subtraction does a good 
job of removing essentially all flux, particularly in K. 

then with the 2MASS-Bessell & Brett filter relation^] we finally 
derive (J - K) , s = 0.70 ± 0.07. 

From our NACO "light curve" using Eq. [2] and our lens 
model , we find after dereddening (J-K)o rS = 0.66 + 0.51. While 
the uncertainty derived from linear regression is large, this inde- 
pendent source color determination is very consistent with the 
colors found by |Gould et al.| ( [20T0| ) as well as those of Bennett 
|et al.| ( [2008] ) and strengthens the case for the source being a K4-5 
dwarf in the Bulge at 7.51 ± 0.25 kpc. However, given the better 
accuracy of the Gould et al. ( 2010| ) source color, we adopt their 
value in the following analysis. 



4.5. Lens/planetary system constraints 

From the MOA-2007-BLG-192 light cur ve the I band sourc e 
flux is well determined to I s = 21.44 + 0.08 ( |Bennett et al.|2008| ). 
Using the source color derived in the previous section and the 
extinction coefficients determined from the IRSF data, we can 
translate this / band estimate into the NACO passbands to derive 
J s = 19.67+0.12,77, = 18.78 + 0.10,/^ = 18.54 + 0.10 (2MASS 
system). 

Using our best lens model and Eq. ^ we then derive the 
following estimates for the apparent lens flux: // = 20.98 + 
0.30, Ht = 19.91 + 0.30,^/ = 19.29 + 0.20 from Epoch 1 and 
// = 20.59 + 0.40, Hi = 20.10 + 0.50, K t = 19.04 + 0.20 from 
Epoch 2. Taking the weighted average we finally get as best esti- 
mate for the lens flux: // = 20.73 + 0.32, Hi = 19.94 + 0.35, K x = 
19.16 + 0.20. 

We now can use mass luminosity relations to translate the 
photometry estimates of the apparent lens flux into est imates of 
the planetary host star mass. We adopt the relations of Del fosse| 
|et al.| 42000] ) for M-dwarfs (with masses > 0.10 M Q ) and 

2 http://www.astro.caltech.edu/ jmc/2mass/v3/transformations/ 



|Chabrier et al.| (|2000) for L-dwarfs (masses < 0.10M o ), where 
the transition between the two relations at ~ 0.10 M Q has been 
linearly interpolated. The best lens model gives an estimate for 
the distance and mass of the lens via the measurement of the par- 
allax tie using Eq. ([I]). In Fig.|7]the implied apparent lens bright- 
ness based on the mass-magnitude relations and our constraint 
on the parallax is plotted as a function of lens mass. All bands 
agree that the lens mass is in the range 0.07 < M^/M Q < 0.10 
with a best estimate of M L /M Q = 0.087 + 0.010, preferring a 
stellar over a sub- stellar host. This is co nsistent with the prev ious 
best estimate of M L /M = 0.06 + 0.04 ( [Bennett et al.|2008| ), but 
which was not able to distinguish between the different host star 
possibilities. This new refined lens mass also affects the inferred 
planetary mass of MOA-2007-BLG-192Lb. This is due to a light 
curve degeneracy between the planetary mass ratio q, and the 
source star radius crossing time t*. The detection of light from 
the lens star means that it must be massive enough to be above 
the hydrogen burning threshold, which constrains t* < 0.05 days 
and rules out the c usp crossing models (models I-P of Table 1 in 
|Bennett et al.| ( |2008"] ) ; the remaining surviving models consistent 
with the NACO data are listed here in the appendix in table |A|l). 
This constraint on t* pushes the mass ratio, q, toward somewhat 
smaller values. As a result the range of allowed planetary masses 
is nearly unchanged. 

The physical parameters of the star-planet system can be esti- 
mated by the same t ype of Markov Chain M onte Carlo (MCMC ) 
calculations used in Bennett et al. ( 2008| ) or Dong et al. ( 2009 ). 
But we now add constraints that the lens star must satisfy the 
JHK mass luminosity relations of Delfos se"et al.| ( |2000| ), un- 
der the assumption that 25% of the dust responsible for the 
extinction of the source star is also in the foreground of the 
lens star plus planet system. The uncertainty in the lens mag- 
nitude is taken to be 0.3 mag in each passband. This accounts 
for the uncertainty in the extinction as well as the uncertainty 
in the Delfosse et al. ( |2000| ) mass-luminosity relations, which 
become large at low masses because of the metallicity depen- 
dence of the minimum stellar mass. The parameter values result- 
ing from this calculation are listed in Table [A] The planet mass 
is now while the host star mass is 0.084 ^q^M© 

and the two dimensional star-planet separation during the event 
is a = 0.66*022 AU. The MCMC lens distance estimate is 

D L = 700^20 P c which 

agrees with our more direct estimate 
of 660^Q°pc. This implies that the lens suffers less than half of 
the total extinction toward the source, and our deri ved lens co lors 
are consistent with a late M spectral type ( Leggett et al.|2010| ) of 
the planetary host. 



4.6. Additional constraints from future high angular 
resolution observations 

Another improvement can be achieved by measuring the ampli- 
tude and direction of the relative proper motion of source and 
lens in combination with the microlensing modeling of the par- 
allax signal caused by the Earth's motion. Such physical mea- 
surements break a model degeneracy in the projected Einstein 
radius f E ( [Bennett et al.|2007||2008) . In the case of MOA-2007- 
BLG-192 the degeneracy is particularly acute because of a gap 
in event coverage, with different equally well-fitting models 
requiring widely different projections and hence directions for 
the relative proper motion, even though the models yield similar 
amplitudes: ~ 5 mas yr" 1 . The measurement of both E and r E 
yields the lens system mass M L = c 2 /(4G) 6 E v E . Ideally, the rel- 
ative lens- source proper motion // re i is measured when detecting 
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Fig. 7. Mass-Magnitude relations for K (top ) , J (middle) and 
H (bottom) bands [CIT system] derived from Delfosse et al. 
( [2000] ) for M-dwarfs (wit h masses > 0.10 Mq) and L-dwarfs 
(masses < 0.10M©) from |Chabrier et al.| ( [2550] ). The transi- 
tion between the two relations at ~ 0.10 M Q has been linearly 
interpolated. The black curves show the most likely range of 
distances for the M OA-2007-BLG-192Lb system as found by 
Bennett et al. (2008) and the horizontal lines marks the estimate 
for the lens flux from the NACO data as well as the upper limit 
of the lens flux from measured lens + source flux for a range of 
f ossible interstellar extinctions. 



both th e lens star and the source star as done by Alcock et al 
( |200TT ). 

The two stars will not be fully resolved for many years. 
However, due to the unique stability of the HST point spread 
function (PSF) it is possible to measure source-lens separations 
(with position angles) much smaller than the width of the PSF. 
This is accomplished by measuring the elongation of the com- 
bined lens-source image due to the fact it is a combination of two 
point source images rather than one. The lens and source stars 
of MOA-2007-BLG-192 will already be about 25 mas apart in 
2012. 

Simulations by Bennett et. al. (2007) show that measure- 
ments of both the amplitude and orientation are possible for 
MOA-2007-BLG-192 already in 2012. These measurements 
combined with our modeling will improve the knowledge of sys- 
tem parameters (masses, orbital separation) to about 10%. The 
key point is that the direction of the elongation will give us a 
measurement of the direction of the relative lens- source proper 
motion yu re i and therefore resolve the remaining parallax degen- 
eracies ( [Ghosh et al.|2004||Bennett et al.|2007||2008] ). 



4. 7. Properties of the planetary system 

The effective temperature of the planet, for the parameters of 
the parent star and orbit separation given above, is 47 for 
an albedo of zero, and 40^K for an albedo of 0.5. Based on 
observations of a tenuous atmosphere (20-60 microbars) of ni- 
trogen on Pluto, the temperature of bright surfaces ices on P luto 
at perihelio n is estimated to be between 35-40 K ( Stansberry & 
|Yelle|1999| ). Thus, if nitrogen were available, the surface of this 
planet might look like that of Pluto on the basis of stellar heating 
alone. However, the large mass of the planet compared with that 
of Pluto necessitates examining the possible role of heat from 
the interior of the planet. The maximum temperature possible 
with zero albedo, 54 K, remains below the pure nitrogen melting 
point of 63 K, and well below the methane melting point of 91 
K. 



The present-day terrestrial heat flow (0.087 W/m 2 ) value is 
about 10 times less than the roughly IW/m 2 deposited by the 
lensing star on its planet at local noon. Thus the average heat 
flow coming from the planet itself will not raise the surface tem- 
perature significantly, even for a fully rocky body three times the 
mass of the Earth (Lunine 2010). Of course, we do not know the 
age of the star; were we to use the Hadean heat flow value for the 
Earth (Hopkins et al.|2008 ) for the 3.2 Earth mass body, the in- 
flux from geothermal heating could exceed the energy received 
from the star. The surface temperature could then be above the 
nitrogen melting point, leading to the possibility of liquid ni- 
trogen lakes or seas if the atmospheric pressure were 0. 1 bar or 
more. The lensing star-planet system is likely older than this, and 
hence the planet's heat flow correspondingly less. 

Because the distribution of heat flow on a terrestrial planet 
can be strongly heterogeneous, one could imagine places on 
the surface with much higher heat flow than the average value 
for the planet, such that temperatures might exceed the melt- 
ing point not just of nitrogen but of methane. Thus, if sufficient 
quantities of these molecules were present, the planet's surface 
might have zones resembling the hydrocarbon lakes and seas of 
Saturn's moon Titan. The possibility of liquid water cannot be 
discounted, but would most likely be below the surface or in 
very restricted, volcanically active, locales. 



D. Kubas, J. P. Beaulieu et al.: NACO constraints MOA-2007-BLG-192Lb 



5. Conclusions 

In this study we have presented the analysis of photometric data 
in the near infrared JHKs bands at two different epochs of plan- 
etary microlensing event MOA-2007-BLG-192 , obtained with 
the AO system NACO mounted o n UT4 at ESQ . Acco rding to 
the best-fit lens models as given in |B ennett et al. | ( [2008] ) the dif- 
ference in the magnification of the source is 0.230 ± 0.015 for 
the two epochs. If the lens is non-luminous this would be the 
expected photometric gradient in our data set in the absence any 
blended light contribution. Our ^-band data, when combined 
with results from previous optical data, are inconsistent with 
such a scenario at 4, 8<x. In fact the data imply that there is a sig- 
nificant amount of blended light at the location of MOA-2007- 
BLG-192. Assuming that this blend is the lens, the data favor a 
scenario in which the lens would be a close-by (660^^°pc) late 
M-dwarf . This is consistent with the estimates for a stellar lens as 
based on constraints from extended source and parallax effects as 
discussed in Bennett et al. (2008). While the data available at the 
time of the discovery paper were consistent with a broad range of 
planetary host masses, the new NACO data presented here sup- 
port the hypothesis of a stellar host for MOA-2007-BLG-192Lb. 
Of course it is conceivable that the detected blend stems not from 
the lens, but either from a stellar companion to the source, the 
lens or an unrelated background star. However the probabilities 
for such scenarios are low and using Ockham's razor the most 
likely explanation is that the lens is an M-dwarf, which implies a 
planetary mass of for MOA-2007-BLG-192Lb, plac- 

ing it among the front row of known least massive cool planets 
in orbit around one of the least massive host stars. 

MOA-2007-BLG-192Lb is a landmark exoplanet discovery 
suggesting that planet formation occurs down to to the very low 
mass end of the stellar pop ulation. This is in ag reement with the 
recent statistical results of |Cassan et al.| (2012) that on average, 
every star in the Milky Way hosts at least one planet. MOA- 
2007-BLG-192 is the first microlensing event for which multi 
epoch AO data has been obtained and demonstrates the use- 
fulness of this technique for microlensing, by constraining the 
physical characteristics of microlensing planetary systems and 
providing important experiences to optimize future AO observa- 
tions, which ideally should be carried out in ToO mode for the 
first epoch, to ensure the source is still significantly magnified. 
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Table A.l. Overview of applied calibrations and transformations between the photometric instrumental systems of IRSF and NACO 
into the 2MASS system. We note that the derived zeropoints for NACO are consistent with zeropoints from NACO based on 
regularly taken standard stars (after correction for atmospheric extinction) and that we do not find a color term between the NACO 
and 2MASS system. 



Band 


photometric calibration 




number of stars used (after last clipping) 


IRSF Single Epoch 


J 


^irsf,2mass = 22.854 + 0.005 + /iRSF,inst - 


- 0.046(/ IRS F,inst " #IRSF,inst) + 0.015 


279 


H 


#IRSF,2MASS - 22.919 + 0.003 + #lRSF,inst 


+ 0.016(/iRSF,inst " #IRSF,inst) + 0.024 


451 


Ks 


^IRSF,2MASS = 22.146 + 0.003 + ^iRSF,inst 


+ 0.010(/ IRS F,inst " ^IRSF,inst) + 0.014 


502 



NACO Epoch 1 



Zeropoints calibrated against IRSF 


using star "1 " (adopted) 


stars "l"+"2" 


J 24.247 + 0.018 + / NAC o,inst 


(24.254 + 0.007) 


H 24.012 ± 0.017 + // N Aco,inst 


(23.987 + 0.028) 


Ks 23.128 + 0.017 + Wo,inst 


(23.105 + 0.020) 



NACO Epoch 2 



Zeropoints calibrated against IRSF 


using star "I " (adopted) 


stars "l"+"2" 


J 24.315 + 0.018 + / NAC o,inst 


(24.345 + 0.030) 


H 24.024 + 0.017 + H NAC0Mt 


(24.030 + 0.060) 


Ks 23.067 + 0.017 + Wo,inst 


(23.091 + 0.023) 



Zeropoints aligned with respect to Epoch I [within 4" of target] 

J 24.289 + 0.019 + / N Aco,inst 8~ 
H 24.036 + 0.012 + # N Aco,inst 10 
Ks 23.1 16 + 0.008 + ff N Aco,inst 10 



Table A.2. Parameter Values and MCMC Uncertainties 



parameter value 2-cr range 

M 0.084 ^^M 0.062-0. 120M o 

m 3.2 + _\\ M e 0.8-14.8M e 

a 0.66^ AU 0.35-3. 17 AU 

D L 0.70^ kpc 0.5-1.4 kpc 

Is 21.59 + 0.07 21.46-21.64 

q l.l + _ l 9 6 Xl0" 4 0.3-5.2 X 10" 4 



Table A.3. This table shows the fit parameters for the 8 distinct planetary models for MOA-2007-BLG-192 consistent with the 
NACO data. t' Q = to - 4240 days, to and uo are the time and distance of the closest approach of the source to the lens center-of-mass. 
q and d are the planet:star mass ratio and separation, and 6 is the angle between the source trajectory and the planet-star axis. I s 
is the best-fit source magnitude, and n E and (p E are the magnitude and angle of the microlensing parallax vector. The units for the 
Einstein radius crossing time, t E , the source radius crossing time, U, and t' Q are days, and all other parameters are dimensionless. 



Fit 


x 1 


t E 




u 


d 


e 


q 






u 


/, 




(f>E 


A 


1121.12 


82.5 


5.442 


0.00309 


0.912 


236.9° 


3.7 x 


10" 


-5 


0.040 


21.61 


1.51 


208.5° 


B 


1121.16 


83.2 


5.442 


0.00306 


1.120 


236.8° 


3.7 x 


10" 


-5 


0.041 


21.62 


1.49 


208.7° 


C 


1121.66 


82.9 


5.456 


-0.00349 


0.807 


105.5° 


3.4 x 


10" 


-4 


0.041 


21.60 


1.47 


209.2° 


D 


1122.08 


83.5 


5.448 


-0.00296 


1.113 


121.9° 


5.8 x 


10" 


-5 


0.041 


21.62 


1.44 


209.8° 


E 


1125.41 


83.2 


5.454 


0.00306 


0.890 


240.1° 


7.6 x 


10" 


-5 


0.048 


21.61 


1.19 


337.5° 


F 


1125.44 


82.5 


5.454 


0.00309 


1.118 


239.9° 


7.5 x 


10" 


-5 


0.049 


21.60 


1.16 


332.3° 


G 


1125.48 


81.8 


5.450 


-0.00313 


0.897 


120.3° 


6.4 x 


10" 


-5 


0.049 


21.60 


1.13 


336.5° 


H 


1125.50 


80.8 


5.450 


-0.00309 


1.110 


120.6° 


6.1 x 


10" 


-5 


0.048 


21.61 


1.19 


337.8° 
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